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Abstract. We study the long-time asymptotics of the probability Pt that 
the Riemann-Liouville fractional Brownian motion with Hurst index H does 

el not escape from a fixed interval [—L,L] up to time t. We show that for any 

' H G]0, 1], for both subdiffusion and superdifTusion regimes, this probability obeys 

In(Pt) ~ —t'^^/L'^, i.e. may decay slower than exponential (subdiffusion) or 
faster than exponential (superdifTusion). This implies that survival probability St 
of particles undergoing fractional Brownian motion in a one-dimensional system 
, with randomly placed traps follows In(S't) ~ — n^/'^t^^/^ as i oo, where n is 

the mean density of traps. 



PACS numbers: 05.40.Ca 



O . 1. Introduction 

Consider stochastic process defined as 

^ ■ * T{H + 1/2) Jo (^-T)i/2-ff' 

' where C(t) is Gaussian, delta-correlated noise and H - the Hurst index, is a real 

number in ]0, 1]. 

' The process in Eq. ([T]) is called the Riemann-Liouville fractional Brownian motion 

OO , (fBm), since Xf^ is the solution of the Langevin equation in which random force 

' is the fractional Riemann-Liouville derivative of Gaussian noise. The fBm was first 

formulated in Kolmogorov's paper and later in the papers of Levy [2 , but first 
systematic analysis belongs to Mandelbrot and Van Ness [3] . 

The fBm X^^ is a continuous-time Gaussian process starting at zero, with mean 



' zero, having the correlation function: 

r^H + 3/2) ^^^U" ' ' 2'ij' 

where s < t and F denotes the Gauss hypergeometric function, and the variance: 

E 



ix."^' 



(3) 



2HT^{H + 1/2)' 
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Equations (O and ^ signify that has self-similar but not independent increments: 
for H > 1/2 the increments of the process are positively correlated and the fBm 
shows a superdiffusive behavior, while for H < 1/2 the increments of the process are 
negatively correlated and one deals with subdiffusion. In case H = 1/2 one recovers 
standard Brownian motion with diffusion coefficient D = 1/2. 

Here we study a simple first-passage problem for the fBm, whose general 
understanding is a basic aspect of stochastic processes [4j [5] . Namely, we analyze the 
asymptotic long-time behavior of the probability Pt that the fBm in Eq. ([T]) does 
not escape from the interval [—L, L] up to time t. We note that contrary to a single 
boundary case, for which several rigorous results are available [6l[7], understanding of 
the two boundaries case is still rather controversial: 

(a) One source of confusion stems from a recent tendency of describing all wealth of 
naturally occurring anomalous diffusive processes in terms of the so-called fractional 
diffusion equation, fractional Fokker-Planck and other fractional differential equations, 
regardless of the origin, intrinsic correlations and physics underlying these processes. 
If non-Markovian fBm-type processes are indeed described by fractional diffusion 
equations |8j, one expects that Pt will have an algebraic tail. This would imply that 
the distribution of the adsorption or first exit time from a fixed interval will not have 
all moments. This is, of course, a rather counterintuitive conclusion. 

(b) Survival of a tagged bead of an infinitely long Rouse polymer chain, (whose 
dynamics is an fBm-type process with H = 1/4), in presence of two adsorbing 
boundaries has been discussed in Ref.[9]. Using a path-integral formulation with an 
exact measure of trajectories of such a bead [10], and adapting a classic method of 
images, it was shown that the bead's survival probability obeys 



i.e., is described by a stretched-exponential function of time. 

(c) Numerical simulations of a tagged particle dynamics in a one-dimensional hard-core 
lattice gas - another fBm-type process with _ff = 1/4, - in presence of two adsorbing 
boundaries [llj . and more recent simulations of dynamics of a tagged bead of a finite 
Rouse chain between two traps [12^, suggested both a faster decay of the survival 
probability: 



We set out to show here that, for < i? < 1, i.e., both in the subdiffusive and 
superdiffusive regimes, the survival probability Pt obeys 



^1/2 

ln(P,)^-^, 



(4) 




(5) 



^2H 

In (Pt) ^ const— -y. 



(6) 



which expression can be rewritten, taking advantage of Eq. ^ , as 



In (Pt) ~ const 



(7) 



Note that our result in Eq.® confirms Eq.(|4]) and contradicts (a) and (c), Eq.®. 
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2. Basic equations 

Since we are concerned with tlie large-t beliavior, it will not matter much how we 
define ({t) - as a continuous in time function or as a discrete process, provided that 
we keep all essential features of noise. We thus divide, at a fixed t, the interval [0, t] 
into N {N ^ 1) small subintervals A, (such that AiV = t), and assume that within the 
fc-th subinterval, k ~ 0,1, . . . , N — I, the noise ((r) is constant and equal to C,k/^fK, 
where {Cfc} are independent random variables with normal distribution N[Q, 1]. 

Then, the Riemann-Liouville fBm in Eq.(IT]) can be written down as a "weighted" 
sum of independent random variables: 



1 



r(i/ + i/2) A^TKAa 

where cr; is a non-random function: 



dT 



N 



(t-r) 



1/2- 



(8) 



1=1 



r(i7 + 3/2) 



lH+l/2 _ _ 



/2 



(9) 



Note that ai is a monotonically decreasing function of / for H < 1/2, a constant for 
H = 1/2, and a monotonically increasing function of / for H > 1/2. 
Now, let recti (X) denote the rectangular function: 







1^1 


<L, 


rectL(X) = < 


1/2, 


X = 


= ±L 




.0, 


1^1 


> L 



(10) 



Representing recti(X) via its Fourier transform: 
dk sin{Lk) 



recti (x) 



■ exp [ikX] , 



(11) 



we may now write down the indicator function llmax \Xl^\ < Lj of the event that 

an A^-step trajectory X/^ did not leave the interval [— i, L] as the following A^-fold 
integral: 



I(max|X/^| <L] = 



N 



dki sin{Lki) 



exp 



' 1=1 



N-l+l 



(12) 



Averaging the latter equation, we have then that the probability P/v of this event is 
given by 



P, 



N 



' — OO >.' — OO 



N 

n 

1=1 



dki sin(Lfci) 



exp 




(13) 
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= h- 


hfc2 + . 


• + few, 


Y2 


= ki- 


h fc2 + . 


• + kN-1 


Y3 


= h- 


hfc2 + . 


■ + kN-2 


Yn 


= h, 







we obtain 



Pn = 



N 



dYi sin 11+1)) 



1=1 



Yi - y,+i 



exp 



1 ^ 



1=1 



(14) 



, Yn+1 = 0. (15) 



Now, it is expedient to use the foUowing integral identity for the sinc-function: 



sin (L {Yi-Yi+^)) ^ 1 
Yi - 11+1 2 



j dXexp[iX(l^ -11+1) 



(16) 



Plugging Eq. (|16|) into Ea. (jl5|) . and performing integrations over {Y/}, we finally arrive 
at the following meaningfuU representation of the survival probability: 



Pn = 



dXi 



L N 



exp 



1=1 



L N 

Y\ dXiexp 

^ 1=1 

N 



2af 

N 



.../ l[dXi exp 

00 J — 00 



-Y,{Xi-Xi^,f /2a\ 
1=1 

N 

-J2{Xi-Xi^,f/2a, 



1=1 



Xo = 0. (17) 



Note that the integrand in Eq. lTTl) is an analog of the Wiener measure for the 
Riemann-Liouville fractional Brownian motion. 



3. Hurst index 1/4 < H < 1/2: fBm as a Brownian motion in an expanding 
cage. 



Change the integration variables Xi = aiXi. Then, Eq. pT]) reads 
/ " ■■•/ 

-Lja\ 'J—L/(T2 L/crjY i—-^ 

where = 0. Next, we represent 



ndxi 
—f= exp 



27r 



1 ^ / ^ 2- 

1 / cr;_i 

9 2^ 



(=1 



N 



1=1 



N 



1=1 



(18) 



(19) 



where 
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F {{xi}) ^2^(1-^) xi^, {xi - + E ( 1 - 



N 



(20) 

1=1 ^ ' i=\ ^ ' 

Now, note that it follows from Eq.® that oi ^ l^^^^^ and (1 — ai-i/ai) ^ l/l as 
I 00. Hence, for any H > and any N, the second sum on the right-hand-side of 
Eq.d^OD is bounded: 



max 



N 



El 



< 



(21) 



On the other hand, an elementary analysis shows that for < H < 1/2 the maximal 
absolute value of the first sum on the right-hand-side of Eq. ([20l) grows with N: 



max 



N 



El 



1=1 



Xi^l {xi - X;_i) 



L^CiN 



1-2H 



(22) 



where Ci is TV-independent constant. 



Consequently, the survival probability P/y obeys the following double-sided 
inequality: 



exp [-2.2 (CiiV^-^H + c)] -^N <Pn < cxp [i^ (CiN^-^" 
in which ^'jv is given explicitly by 

L/cti pL/(J2 pL/a^ j'L/cTN ^ 

— L/(Ji J—L/a^ J—L/<7-i 



c)] ^ 



N, 



N 



ndxi 
—j= exp 

-LjON V 27r 



N . .2 
9 



(23) 



(24) 



with = 0. 

One notices now that Eq. p4|) describes the probability that an iV-step Brownian 
motion trajectory a;;, starting at the origin, does not escape from the interval whose 
boundaries move deterministically away from the origin as ±L/ai, I = l,2,...,iV. 
This classical problem has been extensively studied in the probability theory (see 
Ref.|13j and references therein). A lucid derivation of main results and description of 
different approaches can be also found in Ref. JA\ . 

At sufficiently large N, "^/n obeys [13] (note that we appropriately change the 
notations): 



"^N ~ exp 



2 N 



exp 



exp 



(A7V) 



2H 



16Hr2(ff+l/2) L2 

\ 2" 



a E 



£2 



(25) 



l-2H^ 



Note now that for 1/4 < < 1/2, j.^ in Eq.lgS]) decays faster than exp[-A^ 

and hence, in virtue of the inequality in Eq. p^ . determines the decay of the 
survival probability Pat, which yields the result in Eq.®. 
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4. Hurst index 1/2 < H <1: fBm as a Brownian motion in a shrinking 
cage. 

Consider next the superdifTusive case when 1/2 < H < 1. One notices that here 
ai ~ ^ oo as Z ^ oo, while (1 — ai^i/ai) ~ 1//, and hence, F {{xi}), Eg. ipPl) . 

is bounded by a constant for any N, i.e., 

max|i^ {{xi}) I < 2L'C,, C2 ^ V f 1 - —) <J-_\ (- + —) + 

+ (20) 

Consequently, for 1/2 < < 1, the survival probability obeys: 



exp [-L^C2] ^jv < Pjv < exp [L^Cz] ^tv, (27) 

where vl/^ is defined by Eq. . 

Contrary to the situation discussed in the previous subsection, here, i.e. for 
1/2 < H < 1, describes the probability that an iV-step Brownian motion 
trajectory a;;, commencing at the origin, does not escape from the interval whose 
boundaries move deterministically towards the origin, i.e. that it survives in a 
shrinking cage. 

In this subsection we estimate the long-time asymptotical behavior of 5* at using 
an adiabatic approximation described in Ref. |14j . To ascertain the accuracy of this 
approach, in the next section we will present the results of a more rigorous analysis. 

To define an asymptotic behavior of 5* at, consider the solution of a diffusion 
equation 

9PiX,t) _ 1 d^P{X,t) 

-^^-2A dX^ ' Wi-0)-5(X), (28) 
subject to the boundary conditions 

P{X = ±L{t),t) ^0, (29) 

where L{t) = L/crt and at - VA t"'^^^ /r{H + 1/2), Eq.©. 

The basic idea behind the adiabatic approximation is that, if the cage expands 
or shrinks sufficiently slowly, the density distribution approaches the same form as 
in the fixed-cage case, except that the parameters in this probability distribution 
acquire time dependence to satisfy the moving boundary condition [14j. Within this 
approximation, one takes 



P(X.<)~/«)c„s(^) ./(OcoJf , (30) 



where the amplitude f{t) is to be determined. Substituting Eq. (l5Dl) into Eg. ipS)) . one 
finds 
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df{t) _( 



dt \8AL^{t) 



(31) 



Noticing next that the second term on the right-hand-site of Eq.([3T|) decays faster 



(32) 



than the first 


one, 


we find 








f{t) ~ exp 






/■* dT ' 


~ exp 


r 






8A 


Jo L^{r)\ 


16i7r2(iJH 


-1/2) 



and hence, "^t follows 



*t - / dXP{X,t) = -L(t)exp 



exp 



.2 E 



t 



2H 



16Hr^{H + 1/2) £2 

2" 



(33) 



As in the case 1/4 < H < 1/2, determines the decay of the survival probability 
Pf and hence, we obtain the result in Eq.®. Note also that the exponential term in 
Eq. (|55|) . describing leading behavior in the superdiffusive regime with 1/2 < H < 1, 
is exactly the same as the one in Eg . ((251) . describing subdiffusion, but here it defines 
a faster than exponential decay of the survival probability. 



5. Arbitrary Hurst index, < H < I: Upper and lower bounds. 

In this section, our analysis will be based on the following two keystones: 

/ 1/2 \ H 1 /2 

(i) The probability P I maxjXj ~ | < M) that Brownian motion ~ , starting 



at the origin, does not leave an interval [—M, M] up to time t, t sufficiently large, is 
given by (see, e.g., Ref . [T3] ) : 



P fmax|xf-^/^| < m) - exp 



8M2 



(34) 



(ii) A fundamental property of P/v, Eq. (|17p : Pn — ^W(ci, (72, . . • , (Ttv) is a 
monotonically decreasing function of any variable Ufc, i.e. 

dPN , . 

< 0, < (jfc < oo 35 

dak 

Equation (|35p signifies that replacing any or all au by S(fc), such that Uk < S(fc) for 
any k, we will decrease the survival probability and arrive at the lower bound on Pjv! 
if, on contrary, we will replace one or all Ufc by S(fc), such that ak > S(fc) for any k, 
we will increase the survival probability and obtain an upper bound on Pn- 

We are unaware of any statement similar to (ii) made in the literature, apart of a 
less general " Pascal Principle" p6l . It might be thus instructive to demonstrate first 
its validity. 

Let us first single out terms dependent on Xk and X^-i in Eq. (jl7p writting down 
Pn formally as 
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-L 



N 



n 



dXi 



B{{Xi}) 



X 



X 



dX, 



k 



dXk-i 



B' {Xk-i^Xk^Xk+i) , 



where 



(36) 



AT 



E 

l = l,l^k-l,k,k+l 



and 

B' = cxp 

Note now that, trivially, 



[Xi-Xi^iY 

2af 



{Xi-Xi_^f 



L l=k-l 



2af 



exp 



(Xk-l — Xk_2f' {Xk+l ^ ^k)^ 



<1, 



and hence, Pjv is bounded from above by 
Pn < 



dXi 



n 



Bi{Xi})x 
^ dXk dXk-i 



Xo = 0, 



L V^CTk J-L \/27ra-fe-i 



exp 



{Xk — Xk-i] 



(37) 



(38) 



(39) 



(40) 



Performing integration over Xk and Xk-i, and differentiating both sides of the 
inequahty in Eq. ((40|) with respect to Cfc, we get 



N 



dak 



< 



TTCTfe-l 



1 — exp 



N 

n 

= l,l^k-l,k 



dXi 



2TTai 



^i?({X,})<0(41) 



This proves the inequality in Eq. ([55l 



5.1. Bounds: fBm as a Brownian motion with variance a 



N- 



Suppose that we set in Ea. (|17p all cri equal to ctjv such that the survival probability 
in Eg. pT)) becomes 



P' — 



"A dX, 



exp 



{Xi — Xi-if 



1=1 



2a% 



, ^0 = 0. 



(42) 



In virtue of the fundamental property (ii) of P/v, for subdiffusion, i.e., for H such 
that Q < H < 1/2, we have the following upper bound: 
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Pn < P'^ 



(43) 



since we have replaced ct/, (which is for subdiffusion a monotonically decreasing 
function of by its lowest possible value ctat. On contrary, for superdiffusion we 
have an inverse inequality, i.e., a lower bound, 



Pn > P'f. 



(44) 



since in the superdiffusion case ai is a monotonically increasing function of / and thus 
we have replaced it by its highest possible value. 

(H— 1/2 \ 
max |Xt "'"I < M that 

an iV-step Brownian trajectory, starting at the origin, does not leave an interval 
[— M, M], where M = L/aN- Hence, in virtue of Eq. (|34|l . we have 



P' 



N ~ exp 



'N 



8L' 



N 



exp 



cxp 



(AiV) 



2H 



8r2(iJ+l/2) 



.2 E 



-2H- 



L2 



(45) 



The result in Ea. (|^5|) defines a lower bound on the survival probability Pn in 
the superdiffusive case, and simultaneously, represents an upper bound on P^ for 
subdiffusion, < H < 1/2. 



5.2. Bounds: fBm as a Brownian motion in time t^^ . 

Consider finally a lower and an upper bounds on Pjv which still rely on the fundamental 
property (ii) of the survival probability Pn but also involve a little bit different type 
of arguments. 

In essence, in this subsection we proceed to show that Pn can be bounded by 



/L fL 
Pib{X,t)dX <Pt< / Pub{X,t)dX. 
-L J-L 



(46) 



where Pii,{X,t) and Pub{X,t) obey: 
dPibiX,t) I d^Pib{X,t) 



dt 



2Tib[t) dX^ 



Pib{X,t = 0) = d{X), Tib{t) = 



T^{H + 1/2) 



(47) 



- P„,(;,,, . 0) ^ SiX), T^bit) ^ (48) 



dt 2Tub{t) dX 

which have to be solved subject to the boundary condition 

Pib,ub{X ^±L,t) = Q 



(49) 



In other words, we will show that survival probability Pn can be bounded by survival 
probabilities of Brownian motions evolving in time t^^ . 
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Note, that equations similar to Eas- f^T)) or (pS)) have been akeady proposed in the 
Hterature [T7] as some effective differential equations describing fractional Brownian 
motion. However, despite the fact that Ea. (|T7)) reproduces correctly the variance, 
Eq.Q, the Green's function, and, as we will see, the time-dependence of the survival 
probability, it is not exact and can not reproduce correctly the correlations in the fBm 
process, Eq.(l2]). 

Now, we turn back to our result in Eq. (|17p and notice that in the general case 

< < 1, i.e. for superdiffusion, diffusion and subdiffusion regimes, and for any 

1 > 1, 



,2H 



2Hr^{H +1/2) 



9(1), 



(50) 



where 



gil) ^ I'" - [l - I) 



2H 



(51) 



In fact, the inequality in Eq. ((50|) with g{l) defined by Eq. ljSTj) appears to be a very 
good approximation for ai: while ai and show the same behavior as functions 
of I for sufficiently large I, (Too/S(oo) — 1, they differ only by a few per cent also for 
moderate values of I. 

Now, in virtue of (ii), we have the following lower bound on the survival 
probability Pn' 



Pn > 



M 



L N 

n 



dXi 



(=1 



/2^E(/) 



exp 



L (=1 



2S2(/) 



AI N 

n 



dZ, 



exp 



where Zq = and M = V2HT{H + 1/2)L/A^ 



N 

1=1 



{Zi - Zi^i) 



2{p" -{i-iy") 



(52) 



One notices next that the expression in the second line in Eq. ([52|) defines the 
probability that an iV-step trajectory of Brownian motion, starting at the origin and 
evolving in time T ~ (AN)'^^ does not escape from the interval [— M, A/]. Hence, in 
virtue of Ea. p4p . the survival probability Pn is bounded from below by 



N 



> 



exp 



(AiV) 



2H-\ 



16HT^{H + 1/2) 



exp 



■2 E 



L2 



(53) 



We emphasize that this bound holds for any H g]0, 1] and thus applies to both 
subdiffusion and superdiffusion regimes. Note also that it coincides with the results 
in Eqs.jSSl) and ((33l) . 

Next, note that for any I > 1 and any H e]0, 1], we have 



.2H 



r2(ij + 3/2) 

where g{l) is defined by Eq. ((5T|) . 



5(0, 



(54) 
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Hence, in virtue of (ii), we have the following upper bound on the survival 
probability Pfq: 



Pn < 



L N 
M N 



dXi 



exp 



L (=1 



2S2(/) 



n 



dZ, 



-m'" J-m t\ VMi^" -(i- 1)'^) 



exp 



N 



{Zi - Zi^i) 



^ 2{P" - {I - 1)"^^) 



where Zq = and M = r{H + 3/2)L/A". Consequently, 



(55) 



Pn < exp 



(A7V) 



2H 



8r2(i/ + 3/2) L2 



exp 



2H 



.2 E 



{H + 1/2Y 8 



(56) 



This bound also holds for any H g]0, 1], i.e. for both superdiffusion and subdiffusion. 

Note now that bounds in Eas. ([55)l and appear to be sharper than those 
defined by Eqs. (|^5|) . (l44p and Indeed, for superdiffusion the bound in Eg. ([55)1 is 

higher than the one defined by Eqs. (|44|) and (|45)) . since here 2H > 1. For subdiffusion, 
the upper bound in Eqs. (|^51) and (|l5|) is also worse, i.e. higher, than the one defined 
by Eq.dini) since 2H < 2H/{H + 1/2)^ for H <l/2. 

Therefore, the main result of the present paper can be represented as the following 
double-sided inequality on Pjy: 



2H 



1/2)2 



< - 



L2 



-'e 



In (Pat) < 1, 



(57) 



which holds for any H e]0, 1]. Note that the bounds on the right-hand and on the 
left-hand-side coincide, as they should, for H = 1/2. 



6. Conclusions 



To conclude, we have studied the long-time asymptotical behavior of the probability Pt 
that the Riemann-Liouville fractional Brownian motion with Hurst index H does not 
escape from a fixed interval [~L, L] up to time t. We have shown that Pt obeys 
In(Pi) ^ —t^^ jl?. This result is vahd for any R e]0, 1], for both subdiffusion 
(0 < < 1/2) and superdiffusion (1/2 < < 1) regimes, and consequently, 
the decay may be slower than exponential (subdiffusion) or faster than exponential 
(superdiffusion) . 

This decay law has been obtained by a) showing that for 1/4 < < 1 the 
survival probability Pt of the fBm in presence of two fixed adsorbing boundaries is 
determined by the probability that a Brownian motion does escape from the interval 
with moving boundaries and b) by elaborating upper and lower bounds on Pt which 
show the same time dependence. These bounds stem from some fundamental property 
of the survival probability, Ea. p5p . and controllable approximation of the fractional 
Brownian motion by standard Brownian motion evolving in time T = t^^ . 

The obtained result for the survival probability decay impHes, in particular, that 
the survival probability St of particles undergoing fBm in one-dimensional systems 
with randomly placed traps obeys ln(5't) ^ —r?l^t^^l'^ , where n is the mean density 
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of traps. This expression generalizes famous Balagurov and Vaks [TH] and Donsker 
and Varadhan [12] result over the case of anomalous diffusion described as fractional 
Brownian motion. For H = 1/4, we find the In(S't) '-^ —t^/^ law, which was previously 
obtained in Ref.[n]. 

Finally, we note that the analysis presented in this paper can be straightforwardly 
generalized to fBm taking place in higher-dimensional spaces, the case of Weyl 
fractional Brownian motion, as well as for evaluation of bounds on the distribution 
function of the range of fBm, and on the survival probability of the fBm in presence 
of one-sided or two-sided moving boundaries. 
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